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ABSTRACT 

We  present  an  analysis  of  the  Far  Ultraviolet  Spectroscopic  Explorer  (FUSE)  spectra  of  the  little-known  southern 
nova-like  cataclysmic  variable,  BB  Doradus.  The  spectrum  was  obtained  as  part  of  our  Cycle  8  FUSE  survey  of  high- 
declination  nova-like  stars.  The  FUSE  spectrum  of  BB  Dor,  observed  in  a  high  state,  is  modeled  with  an  accretion 
disk  with  a  very  low  inclination  (possibly  lower  than  10°).  Assuming  an  average  white  dwarf  (WD)  mass  of  0.8  Mq 
leads  to  a  mass  accretion  rate  of  10^®  Mq  yr^'  and  a  distance  on  the  order  of  ~650  pc,  consistent  with  the  extremely 
low  Galactic  reddening  in  the  direction  of  BB  Dor.  The  spectrum  presents  some  broad  and  deep  silicon  and  sulfur  ab¬ 
sorption  lines,  indicating  that  these  elements  are  overabundant  by  3  and  20  times  solar,  respectively. 

Subject  headings:  accretion,  accretion  disks  —  novae,  cataclysmic  variables  —  stars:  abundances  — 
stars:  individual  (BB  Doradus)  —  white  dwarfs 

Online  material:  color  figures 


1.  INTRODUCTION 
1.1.  Nova-like  Cataclysmic  Variables 

Cataclysmic  variables  (CVs)  are  close  binaries  in  which  the 
primary,  a  white  dwarf  ( WD),  accretes  matter  and  angular  mo¬ 
mentum  from  the  secondary,  a  main-sequence  star,  filling  its  Roche 
lobe.  The  matter  is  transferred  by  means  of  either  an  accretion 
disk  around  the  WD  or  an  accretion  column,  formed  when  the 
WD  has  a  strong  (~10  MG)  magnetic  field.  Ongoing  accretion 
at  a  low  rate  (quiescence)  is  interrupted  every  few  weeks  to  months 
by  intense  accretion  (outburst)  of  days  to  weeks  (a  dwarf  nova 
accretion  event),  and  every  few  thousand  years  by  a  thermonuclear 
explosion  (TNR;  the  classical  nova  event).  CV  systems  are  divided 
into  subclasses  according  to  the  duration,  occurrence,  and  ampli¬ 
tude  of  their  outbursts:  dwarfnovae  (DNs)  are  found  mostly  in  the 
low  state,  while  nova-like  CVs  (NLs)  exhibit  the  spectroscopic 
and  photometric  characteristics  of  novae  between  or  after  out¬ 
bursts,  but  have  never  had  a  recorded  outburst.  NLs  form  a  less 
homogeneous  class  and  are  divided  into  subtypes  depending  on 
certain  properties,  such  as  eclipsing  systems,  magnetic  systems, 
or  systems  that  go  into  unexpected  low  states  (Warner  1995). 
Nonmagnetic  NLs  can  be  divided  into  UX  UMa  systems,  which 
remain  in  a  state  of  high  optical  brightness  or  a  “permanent  out- 

’  Based  on  observations  made  with  the  NASA-CNES-CSA  Far  Ultraviolet 
Spectroscopic  Explorer,  which  is  operated  for  NASA  by  the  Johns  Hopkins 
University  under  NASA  contract  NAS5-32985. 
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burst  state,”  and  VY  Scl  systems  (anti-dwarf  novae),  which 
experience  unexpected  low  states  when  the  optical  brightness 
plummets. 

CVs  are  believed  to  evolve  from  long  period  (~2  days)  to¬ 
ward  short  period  (^  1  hr)  as  the  mass  ratio  q  decreases,  crossing 
the  2-3  hr  “period  gap,”  where  few  systems  are  found.  Theory 
predicts  that  the  period  of  those  systems  with  a  sufficiently  low  q 
will  increase  as  q  continues  to  decrease,  producing  the  so-called 
“period  bouncers”  (Patterson  et  al.  2005a).  The  driving  mech¬ 
anism  behind  the  mass  transfer  (M)  and  hence,  the  evolution  of 
CVs,  is  believed  to  be  angular  momentum  loss,  dominated  by 
magnetic  stellar  wind  braking  at  periods  above  3  hr  and  by  grav¬ 
itational  wave  emission  below  this  period.  The  secondaries  are 
stripped  to  <0.08  Mq  on  a  timescale  of  only  IM  Gyr  after  they 
form  as  close  binaries  above  the  gap.  Consequently,  in  the  life¬ 
time  of  the  Galaxy,  the  vast  majority  of  CVs  should  have  evolved 
to  a  period  minimum  near  80  minutes,  and  should  now  have  de¬ 
generate  brown  dwarflike  secondaries  (Howell  et  al.  2001).  How¬ 
ever,  recent  models  (Andronov  et  al.  2003)  suggest  (1)  much 
lower  angular  momentum  loss  rates,  such  that  it  takes  10-12  Gyr 
to  even  reach  a  3  hr  period,  and  (2)  multiple  evolutionary  tracks 
yielding  different  populations  of  CVs  above  and  below  the  3  hr 
orbital  period. 

The  accretion  disks  and  underlying  accreting  WDs  in  CVs  can 
provide  crucial  clues  and  constraints  on  these  evolutionary  sce¬ 
narios.  It  is  important  to  derive  M  from  observations,  as  the  mass 
transfer  rate  is  tied  to  the  angular  momentum  loss.  As  to  the  WDs, 
they  are  the  central  engines  of  the  observed  outbursts,  either  as  po¬ 
tential  wells  for  the  release  of  accretional  energy  or  as  the  sites  of 
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explosive  TNR  shell  burning.  In  addition,  the  differences  be¬ 
tween  the  temperatures,  rotation,  and  chemical  abundances  of 
CV  WDs  and  single  isolated  WDs  provide  clues  as  to  the  ef¬ 
fects  of  accretion,  diffusion,  and  long-term  heating  and  evolu¬ 
tion  (Sion  1995). 

With  a  wavelength  range  of  905-1187  A,  the  Far  Ultraviolet 
Spectroscopic  Explorer  {FUSE)  covers  that  part  of  the  far-ultraviolet 
(FUV)  spectram  where  the  hotter  CV  component  is  dominant. 
For  NLs,  the  dominant  component  in  the  FUSE  range  is  usually 
the  accretion  disk.  Additional  possible  components  are  the  hot 
white  dwarf,  the  boundary  layer  or  accretion  belt,  and  other  hot 
regions  on  or  close  to  the  WD.  Many  CV  FUV  spectra  are  mod¬ 
eled  with  two  components,  usually  a  WD  plus  a  disk.  For  each 
system,  one  usually  derives  from  the  FUSE  spectram  the  temper¬ 
ature  of  the  WD  (Teff),  the  gravity,  the  rotational  velocity  (Uot  sin  i), 
the  chemical  abundances,  the  mass  accretion  rate,  the  inclination, 
and  the  distance  to  the  system. 

It  has  therefore  been  very  important  to  try  to  obtain  FUSE 
data  for  more  CV  systems,  in  order  to  populate  the  Teff-period 
and  the  A/-period  parameter  spaces.  There  is  a  need  to  enlarge 
the  sample  to  obtain  an  accurate  global  picture  of  the  systems 
above  and  below  the  period  gap.  For  that  reason,  we  proposed 
to  observe  with  FUSE  a  set  of  1 6  high-declination  dwarf  novae 
(our  Cycle  7  survey  of  DNs)  and  16  high-declination  NLs  (our 
Cycle  8  survey  of  NLs,  all  chosen  from  the  online  CV  catalog  of 
Downes  et  ah).  Unfortunately,  our  FUSE  NL  survey  was  cut 
short  due  to  the  fatal  failure  of  the  reaction  wheel  of  the  tele¬ 
scope.  Of  the  16  targets,  only  two  were  observed:  BB  Dor  and 
P83 1-57.  The  analysis  of  the  FUSE  spectram  of  P83 1-57  is  pre¬ 
sented  elsewhere  (P.  Barrett  et  al.  2008,  in  preparation).  In  this 
work,  we  present  the  result  of  the  FUSE  spectroscopic  analysis 
of  BB  Dor. 

1.2.  The  Nova-like  CV  BB  Doradus 

BB  Dor  (also  known  as  EC  05287—5847,  an  object  from  the 
Edinburgh-Cape  Survey  [Chen  et  al.  2001])  was  spectroscop¬ 
ically  identified  as  a  CV  in  1987  December,  and  was  seen  fainter 
than  V  ^  16.5  until  1992  November  (Chen  et  al.  2001).  Since 
that  time,  it  seems  to  be  in  bright  state  {V  ^  14.6-13.6).  A  first 
(extremely  uncertain)  estimate  of  its  period  by  Chen  et  al.  (2001) 
put  it  right  in  the  middle  of  the  period  gap,  with  a  period  of 
0.107  days,  or  2.57  hr.  However,  a  45  day  observation  (Patterson 
2002;  Patterson  et  al.  2005b)  corrected  this  value  to  0. 14923  days, 
or  3.58  hr,  as  expected  for  VY  Scl  stars,  which  usually  have  pe¬ 
riods  between  ~3  and  ^4  hr.  BB  Dor  exhibits  superhumps,  with 
a  period  excess  of  £  =  9.4%  (Patterson  2002);  it  has  the  largest 
known  fractional  superhump  of  all  systems,  and  therefore  pro¬ 
vides  a  calibrating  point  at  large  e  (Patterson  et  al.  2005b). 
BB  Dor  is  a  VY  Scl  NL,  and  so  far  no  magnetic  activity  has  been 
found.  The  system  also  exhibits  quasi-periodic  oscillations  (Chen 
et  al.  2001).  Chen  et  al.  (2001)  also  noted  that  due  to  its  very 
narrow  Balmer  lines  (with  wings  extending  ^650  km  s^^  from 
the  line  center),  BB  Dor  is  most  likely  a  low-inclination  system. 
From  the  Austral  Variable  Star  Observer  Network  (AVSON),  it 
appears  that  BB  Dor  has  been  around  V  «  13.6-14.6  mag  (with 
an  error  <0. 1  mag)  for  the  last  couple  of  years  and  can  vary  as 
much  as  0.4  mag  within  1  hr.  The  reddening  value  E(B  —  V) 
toward  BB  Dor  is  not  known.  However,  the  Galactic  reddening 
value  in  the  direction  of  BB  Dor  as  inferred  from  the  100  pm  dust 
emission  map  (Schlegel  et  al.  1998)  is  very  small:  E{B  —  V)  = 
~  0.03.  Since  the  dust  map  gives  the  maximum  dust  extinc¬ 
tion,  we  deduce  that  BB  Dor  must  have  a  negligible  reddening 
[E(B  —  V)  0.01  almost],  which  should  not  affect  the  results; 
consequently  we  do  not  deredden  its  FUSE  spectrum. 


2.  OBSERVATIONS 
2.1.  The  FUSE  Spectrum  of  BB  Dor 

BB  Doradus  was  observed  with  FUSE  on  2007  July  9  at 
02:51:09  (UT),  and  three  exposures  were  obtained  (H9030301; 
exposures  001,  002,  003).  The  FUSE  spectrum  was  obtained 
through  the  low-resolution  (LWRS;  30"  x  30")  aperture  in  TIME 
TAG  mode.  The  data  were  processed  with  the  last  and  final  ver¬ 
sion  of  CalFUSE,  version  3.2.0  (Dixon  et  al.  2007). 

A  look  at  the  count  rate  graph  (generated  by  CalFUSE)  in¬ 
dicates  that  the  data  for  each  exposure  were  collected  mainly 
during  the  first  1000  s,  followed  by  a  jitter  lasting  about  1500  s. 
Toward  the  end  of  each  exposure,  the  count  rate  dropped  to  zero. 
Therefore,  while  the  coadded  exposures  for  the  8  individual  seg¬ 
ments  after  processing  by  CalFUSE  give  a  good  exposure  time 
(after  screening)  of  between  about  8750  s  for  segments  SiCla, 
SiClb,  LiFla,  and  LiFlb  (segments  1)  and  9750  s  for  segments 
SiC2a,  SiC2b,  LiF2a,  and  LiF2b  (segments  2),  the  count  rate  plot 
indicates  that  the  total  good  collection  time  ranges  between  3350  s 
(segments  1)  and  3650  s  (segments  2).  We  therefore  weighted  these 
segments  accordingly,  multiplying  the  flux  (on  the  order  of  1  x 
10^*^  erg  s^'cm  in  the  online  preview)  by  a  factor  of 

2.61  (segments  1)  and  2.67  (segments  2).  After  that,  we  pro¬ 
cessed  the  data  as  usual  (see  below). 

The  spectral  regions  covered  by  the  spectral  channels  overlap, 
and  these  overlap  regions  are  then  used  to  scale  the  spectra  in  the 
SiCl,  SiC2,  and  LiF2  channels  to  the  flux  in  the  LiFl  channel. 
The  low-sensitivity  portions  (usually  the  edges)  of  each  channel 
are  discarded.  In  the  present  case,  the  SiC2b  spectral  segment 
was  very  noisy,  and  we  discarded  it  too.  We  also  carried  out  a 
visual  inspection  of  the  FUSE  channels  to  locate  “the  worm” 
(a  wire  flaw  in  the  FUSE  spectrograph  that  occults  the  saggital 
focus),  and  we  manually  discarded  those  portions  of  the  spec¬ 
tram  affected  by  the  worm.  We  combined  the  individual  expo¬ 
sures  and  channels  to  create  a  time-averaged  spectram  weighting 
the  flux  in  each  output  datum  by  the  exposure  time  and  sensitiv¬ 
ity  of  the  input  exposure  and  channel  of  origin.  The  final  product 
is  a  spectram  that  covers  the  full  FUSE  wavelength  range  of 
905-1 1 87  A.  Because  we  disregarded  the  right  edge  of  both  the 
Sic  1  a  and  SiC2b  segments,  there  is  a  gap  between  the  LiF  1  a  and 
LiF2b  segments  (~  1082-^  1087  A). 

The  FUSE  spectram  of  BB  Dor  is  shown  in  Figure  1 .  Note  that 
the  FUSE  exposures  of  BB  Dor  consist  entirely  of  observations 
carried  out  during  the  day.  The  short  night  exposure  was  ignored 
because  of  jitters  and  other  problems.  We  suspect  that  most  of 
the  emission  lines  are  due  to  helio-  and  geocoronal  emissions. 

2.2.  The  FUSE  Lines 

The  FUSE  spectra  of  CVs  exhibit  mainly  broad  absorption 
lines  from  the  accretion  disk  and  WD,  as  well  as  sharp  absorption 
lines  from  circumstellar  (or  circumbinary)  material  and  the  inter¬ 
stellar  medium  (ISM).  Emission  lines  from  the  hotter  regions 
are  usually  broadened,  depending  on  the  inclination  angle  i,  due 
to  the  large  Keplerian  velocity  in  the  inner  disk.  Sharp  emission 
lines  from  air  glow  (geo-  and  heliocoronal  in  origin)  are  also 
present,  due  mainly  to  the  reverberation  of  sunlight  inside  the 
FUSE  telescope  during  daytime  observations.  The  present  F’t/SZ? 
spectram  of  BB  Dor  presents  such  a  complexity,  and  all  the  lines 
are  listed  in  Table  1 . 

Absorption  lines. —  The  main  characteristic  of  the  FUSE 
spectram  of  BB  Dor  is  the  broad  Ly/3  absorption  feature  due  to 
either  the  exposed  WD,  the  disk  at  low  inclination  angle  i,  or 
possibly  both  (see  §  4).  Absorption  features  from  the  higher 
orders  of  the  Lyman  series  are  also  clearly  visible,  indicating  a 
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Fig.  1 . —  Solar  abundance,  single-WD  model.  The  best-fit  WD  model  assuming  solar  abundances  (solid  black  line)  is  shown  together  with  the  FUSE  spectrum  of 
BB  Dor  (red  line).  Blue  lines  show  the  segments  of  the  spectrum  that  were  masked  before  the  fitting.  The  WD  has  a  temperature  T  =  37,000  K,  a  projected  rotational 
velocity  Trot  sin  i  =  400  km  s“*,  and  log  g  =  8.3.  The  distance  obtained  is  rf  =  217  pc,  and  =  0.3348. 


temperature  T  ^  35, 000  K.  The  other  main  absorption  features 
expected  at  this  temperature  are  C  iii  (1175  A),  Si  iii  (fallOS- 
1114  A  and  sal  140-1 144  A),  Si  iv  («1 120-1 130  A),  and  C  ii 
(1010  A).  The  spectrum  is,  however,  also  marked  by  some  high- 
order  ionization  absorption  lines,  such  as  S  iv  (1063  and  1073  A), 
Si  IV  (1066.6  A),  and  the  O  vi  doublet.  The  S  iv  (1006  A),  C  ii 
(1010  A),  andS  IV  ( 1 099  A)  absorption  multiplets  are  all  unre¬ 
solved  and  shallow.  However,  we  have  indicated  them  on  the 
figures,  as  they  appear  deeper  in  the  modeling.  On  the  other  hand, 
the  S  IV  (1063  and  1073  A)  and  Si  iii  (^1010  A)  absorption  lines 
appear  deeper  in  the  observed  spectmm  than  in  the  modeling.  The 
presence  of  high-order  ionization  species  (such  as  the  oxygen 
doublet)  indicates  the  presence  of  a  hotter  absorbing  component 
above  the  main  FUV-emitting  region.  The  P  ii  line  (~961  A)  has 
been  marked  in  the  spectrum,  although  the  feature  is  most  prob¬ 
ably  due  to  noise.  All  the  lines  are  listed  in  Table  1  with  their 
wavelengths.  Most  of  the  broad  absorption  lines  and  features 
associated  with  the  source  are  redshifted  by  about  0.6-0. 9  A, 
corresponding  to  a  receding  velocity  of  200  ±  40  km  s^  * .  It  seems 
very  likely  that  this  redshift  is  due  to  the  orbital  motion  of  the 
WD;  however,  we  do  not  have  enough  time  tag  data  to  verify  this 
assumption. 

Emission  lines. —  The  spectrum  of  BB  Dor  exhibits  sharp 
emission  lines,  including  all  the  orders  of  the  Lyman  series,  N  iv 
(-923  A),  S  VI  (933.5  and  944.5  A),  C  iii  (977  A),  He  ii  (992  A), 
and  He  I  (1 068  A).  Some  of  these  lines,  in  particular  the  H  i  lines,  the 
He  I  and  He  ii  lines,  the  C  iii  (977  A)  line,  and  the  O  vi  doublet, 
seem  to  be  due  to  heliocoronal  emission  (sunlight  reflected  in¬ 
side  the  telescope),  which  contaminates  the  SiC  channels.  The 
segments  of  the  spectrum  from  the  LiF  channels  overlapping 


the  Sic  Channels  down  to  about  1 000  A  do  not  show  any  oxygen 
lines,  and  we  therefore  did  not  include  these  spectral  regions  of 
the  Sic  channels.  However,  other  emission  lines  may  be  solely 
from  the  source.  In  many  NLs,  the  emission  lines  from  the  disk 
are  broadened  by  the  Keplerian  velocity  and  are  consequently 
easily  identified.  However,  Chen  et  al.  (2001)  have  detected  very 
narrow  Balmer  lines  and  concluded  that  BB  Dor  is  most  likely  a 
low-inclination  system.  Because  of  this,  the  identification  of  the 
emission  lines  from  the  source  is  not  trivial.  This  issue  is  further 
complicated  by  the  lack  of  night  exposure  and  the  very  low  signal- 
to-noise  ratio,  which  makes  it  dilficult  to  analyze  the  two- 
dimensional  image  of  the  spectrum  for  the  effect  of  scattered 
sunlight  on  the  emission  lines. 

Since  the  system  is  a  NL  in  a  high  state,  we  first  expect  the  N  iv 
(«923  A),  S  VI  (933.5  and  944.5  A),  and  C  iii  (977  A)  emission 
lines  to  be  from  the  source,  with  some  possible  heliocoronal  con¬ 
tamination  mostly  aifecting  the  carbon  emission  lines.  The  N  iv 
(923.06  A)  emission  line  is  contaminated  with  the  H  i  (923. 1 5  A) 
line,  while  the  N  iv  (922.52  and  924.28  A)  lines  are  not.  How¬ 
ever,  the  complete  absence  of  narrow  emission  lines  from  the  oxy¬ 
gen  doublet  and  C  iii  (1175  A)  seems  to  indicate  otherwise;  i.e., 
all  the  sharp  emission  lines  might  be  heliocoronal  in  origin.  To 
further  check  this  possibility,  we  measured  the  relative  intensities 
of  the  narrow  emission  lines  below  1000  A  and  compared  them 
with  those  of  the  solar  spectrum  (Curdt  et  al.  2001).  We  found 
good  agreement  with  the  solar  disk-average  quiet-Sun  data,  ex¬ 
cept  for  S  VI  (944. 5  A),  which  has  a  higher  relative  intensity  in  the 
spectrum  of  BB  Dor.  This  could  be  due  to  sunspot  activity  during 
the  FUSE  observation  of  BB  Dor,  as  this  line  intensity  increases 
by  a  factor  of  30  inside  sunspots,  while  the  other  lines  increase 
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TABLE  1 

FUSE  Line  Identification 


Ion 

'^rest 

(A) 

■^‘•obs 

(A) 

Origin/Nature  of  Line® 

Hi . 

918.13 

918.0 

c/e 

919.35 

919.3 

c/e 

920.96 

920.9 

c/e 

N  IV . 

922.52 

922.3 

s/e 

923.06 

923.0 

s,  c/e 

Hi . 

923.15 

923.2 

c/e 

N  IV . 

923.22 

923.2 

c/e 

924.28 

924.2 

c/e 

Hi . 

926.23 

926.2 

c/e 

930.75 

930.7 

c/e 

S  VI . 

933.50 

933.4 

c/e 

Hi . 

937.80 

937.8 

c/e 

S  VI . 

944.50 

944.4 

c/e 

Hi . 

949.74 

949.7 

c/e 

Ni . 

952.40 

925.5 

ism /a 

953.42 

953.66 

953.97 

954.10 

954.2 

ism /a 

P  ii . 

961.04 

961.2 

ism /a? 

Hi . 

972.54 

972.5 

c/e 

Cm . 

977.02 

977.0 

c/e 

Si  III . 

993.52 

993.2 

s?/a 

994.19 

994.9 

s?/a 

997.39 

997.1 

s?/a 

S  IV . 

1006.07 

s/a  -  unresolved 

1006.39 

s/a  -  unresolved 

C  II . 

1009.85 

s/a  -  unresolved 

1010.08 

s/a  -  unresolved 

1010.37 

s/a  -  unresolved 

Hi . 

1025.77 

1025.7 

c/e 

0  VI . 

1031.91 

1032.9 

s/a 

C  II . 

1036.34 

1036.4 

ism /a 

0  VI . 

1037.61 

1038.3 

s/a 

0  I . 

1039.10 

1039.3 

ism /a 

An . 

1048.20 

1048.3 

ism /a 

S  IV . 

1062.65 

1063.4 

s/a 

Si  IV . 

1066.60 

1067.1 

s/a 

An . 

1066.66 

1066.8 

ism /a 

S  IV . 

1072.97 

1073.9 

s/a 

1073.52 

1073.9 

s/a 

S  IV . 

1098.36 

s/a  -  unresolved 

1098.93 

s/a  -  unresolved 

1099.48 

s/a  -  unresolved 

1100.53 

s/a  -  unresolved 

Si  III . 

1108.36 

1108.8 

s/a 

1109.94 

1110.6 

s/a 

~1113.20 

1113.8 

s/a 

Si  IV . 

1122.49 

1123.0 

s/a 

~1128.33 

1128.9 

s/a 

Ni . 

1134.16 

1134.3 

c/a 

1134.42 

1134.3 

c/a 

1134.98 

1135.2 

c/a 

S  i . 

1145.10 

1145.0 

ism /a 

He  II . 

1168.61 

1168.6 

c/e 

Cm . 

1174.90 

1176.4 

s/a  -  unresolved 

1175.26 

1176.4 

s/a  -  unresolved 

1175.60 

1176.4 

s/a  -  unresolved 

1175.71 

1176.4 

s/a  -  unresolved 

1176.00 

1176.4 

s/a  -  unresolved 

1176.40 

1176.4 

s/a  -  unresolved 

^  (a):  absorption;  (e): 

emission;  (c):  contamination  (e.g.,  air  glow;  heliocoronal 

emission);  (s):  source;  ‘ 

ism”indicates  interstellar  medium. 

by  at  most  a  factor  of  3 .  We  therefore  conclude  that  all  the  sharp 
emission  lines  in  the  FUSE  spectrum  of  BB  Dor  are  due  to  sun¬ 
light  reflected  inside  the  telescope. 

3.  SPECTRAL  MODELING 

We  created  a  grid  of  models  of  synthetic  spectra  of  WDs  and 
accretion  disks  for  different  values  of  the  WD  temperature 
gravity  log  g,  projected  rotational  velocity  V^ot  sin  i,  inclination  i, 
mass  accretion  rate  M,  and  abundances.  We  then  ran  a  fit" 
ting  program  to  And  the  best  At  for  (1)  a  single  WD  component, 
(2)  a  single  accretion  disk  component,  and  (3)  a  combined 
WD-i-accretion  disk.  We  describe  below  how  we  generate  these 
synthetic  spectra  and  how  we  perform  the  fitting. 

We  create  the  synthetic  model  spectra  for  high-gravity  stellar 
atmospheres  using  the  TLUSTY  and  SYNSPEC^  codes  (Hubeny 
1988;  Hubeny  &  Lanz  1995).  Atmospheric  structure  is  com¬ 
puted  (using  TLUSTY)  assuming  a  H-He  LTE  atmosphere;  the 
other  species  are  then  added  in  the  spectrum  synthesis  stage 
using  SYNSPEC.  We  generate  photospheric  models  with  effec¬ 
tive  temperatures  ranging  from  ~20,000-50,000  K  in  increments 
of  1 000  K.  We  chose  values  of  log  g  ranging  between  7.5  and  9.0. 
We  also  varied  the  stellar  rotational  velocity  V^ot  sin;  irom  100 
to  1000  km  s^^  in  steps  of  100  km  s^*  (or  smaller  if  needed). 
The  WD  rotation  ( Vrot  sin ;)  rate  is  determined  by  fitting  the 
WD  model  to  the  spectrum  while  paying  careful  attention  to  the 
line  profiles  in  the  FUSE  spectrum.  We  do  not  carry  out  separate 
fits  to  individual  lines,  but  rather  try  to  fit  the  lines  and  continuum 
simultaneously,  while  paying  careful  attention  to  the  absorption 
lines.  It  is  important  to  note  that  the  depth  of  the  absorption  fea¬ 
tures  depends  not  only  on  the  abundances,  but  also  on  the  ro¬ 
tational  velocity.  Increasing  the  rotational  velocity  reduces  the 
central  depths  of  the  absorption  features,  thus  reducing  the  abun¬ 
dance.  However,  the  widths  of  the  absorption  features  also  in¬ 
crease  with  increasing  rotational  velocity.  As  a  consequence, 
abundances  and  rotational  velocity  effects  are  intertwined,  and 
cannot  always  be  easily  separated  without  additional  informa¬ 
tion  about  one  of  the  parameters.  For  any  WD  mass,  there  is  a 
corresponding  radius,  or  equivalently  one  single  value  of  log  g 
(e.g.,  see  the  mass  radius  relation  from  Hamada  &  Salpeter 
[1961],  or  see  Wood  [1990]  or  Panel  et  al.  [2000]  for  WDs  of 
different  compositions  and  nonzero  temperature).  Therefore,  by 
scaling  the  theoretical  spectrum  to  the  observed  one,  we  obtain 
the  distance  to  the  system. 

We  model  accretion  disk  spectra  by  first  assuming  that  the 
disk  is  made  of  a  collection  of  annuli,  where  each  annulus  has  a 
temperature  T{r)  and  gravity  log  g{r),  as  given  by  the  standard 
disk  model  (Shakura  &  Sunyaev  1973;  Pringle  1981),  for  a 
given  central  mass  AAvd  and  accretion  rate  M.  A  variant  of  the 
TLUSTY  code,  called  TLDISK,  is  then  used,  which  generates 
an  atmosphere  model  for  each  annulus  that  is  then  used  as  input 
for  SYNSPEC.  The  contribution  of  all  the  annuli  are  then  com¬ 
bined  using  DISKSYN,  and  a  final  spectrum  is  obtained  for  any 
given  inclination  angle.  A  detailed  explanation  of  the  procedure 
is  given  in  Wade  &  Hubeny  (1998).  In  the  present  work,  we  do 
not  use  the  grid  of  synthetic  accretion  disk  spectra  tabulated  by 
Wade  &  Hubeny  (1998);  instead,  we  generate  them.  This  allows 
us  to  compute  disk  spectra  assuming  nonsolar  abundances  and 
for  any  inclination  angle  (the  disk  spectra  of  Wade  &  Hubeny 
[1998]  have  solar  abundances  and  have  been  generated  for  a 
specific  value  of  the  inclination  ;). 

Before  carrying  out  a  synthetic  spectral  fit  of  the  spectra,  we 
masked  portions  of  the  spectra  with  strong  emission  lines,  strong 

^  TLUSTY,  ver.  200  and  SYNSPEC,  ver.  48;  see  http://nova.astro.umd.edu. 
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Fig.  2. — Nonsolar  abundance,  single-WD  model.  The  best-fit  WD  model  (solid  black  line)  to  the  FUSE  spectmm  of  BB  Dor  (solid  gray  line)  is  shown,  where  the  abun¬ 
dances  of  sulfur  and  silicon  have  been  set  to  [S]  =  20  solar  and  [Si]  =  3  solar.  The  abundances  of  all  the  other  elements  have  been  set  to  solar.  Dashed  gray  line  segments 
show  the  segments  of  the  spectrum  that  were  masked  before  the  fitting.  The  WD  has  a  temperature  T  =  37,000  K,  a  projected  rotational  velocity  sin  i  =  400  km  s“*, 
and  logg  =  8.3.  The  distance  obtained  is  =  21 1  pc,  and  =  0.3129.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 


ISM  absorption  lines,  detector  noise,  and  air  glow.  The  regions 
excluded  from  the  fit  appear  in  blue  in  Figure  1. 

After  generating  grids  of  models  for  the  FUSE  spectrum  of 
BB  Dor,  we  use  FIT  (Press  et  al.  1992),  a  minimization  rou¬ 
tine,  to  compute  the  reduced  (i-e.,  x^  per  number  of  degrees 
of  freedom)  and  scale  factor  (which  gives  the  distance)  for  each 
model  fit.  While  we  use  a  x^  minimization  technique,  we  do  not 
blindly  select  the  least-x^  models,  but  we  also  examine  the  mod¬ 
els  that  best  fit  some  of  the  features  such  as  absorption  lines  and, 
when  possible,  the  slope  of  the  wings  of  the  broad  Lyman  absorp¬ 
tion  features. 

Initially,  we  generate  solar  abundance  models,  and  when  a 
good  fit  is  found,  we  start  varying  the  chemical  abundances  of 
C,  N,  S,  and  Si  to  fit  the  absorption  features  of  the  spectrum.  In 
particular,  the  carbon  abundance  was  set  using  the  C  ii  (1010  A) 
and  C  III  ( 1 1 75  A)  multiplets,  the  sulfur  abundance  was  set  using 
the  S IV  (1063  and  1073  A)  lines,  and  the  silicon  abundance  was 
set  using  the  Si  iv  (1067 , 1023  ,  and  1028  A)  and  Si  iii  (~1110 
and«1138-1146  A)  lines. 

4.  RESULTS 

The  data  obtained  by  AVSON  imply  that  BB  Dor  has  been  in 
a  high  state  around  V  «  13.6-14.6  for  the  last  couple  of  years. 
We  therefore  expect  the  FUSE  spectrum  to  be  dominated  by 
flux  from  the  accretion  disk  with  a  relatively  high  M.  However, 
in  our  modeling  we  follow  a  systematic  approach  that  consists 
of  fitting  (1)  a  single  WD,  (2)  a  single  accretion  disk,  and  (3)  a 
WD+accretion  disk  composite. 


4. 1 .  White  Dwarf 

Since  we  do  not  have  any  information  about  the  mass  of  the 
WD,  we  look  for  all  the  best-fit  models  in  the  log  g-T^s  plane. 
Namely,  for  each  assumed  value  of  log  g,  we  vary  the  tempera¬ 
ture  to  find  the  best-fit  model.  As  expected,  we  find  that  the  tem¬ 
perature  is  somewhere  between  34,000  K  (for  log  g  =  7.5)  and 
40,000  K  (for  log  g  =  9.0),  with  a  distance  between  377  and 
135  pc,  respectively  (see  Table  1).  The  least  xl  is  obtained  for 
log  g  =  8.0-8.65.  We  chose  the  intermediate-value  log  <7  =  8.3 
model,  with  T  =  37,000  K,  J  =  247  pc,  and  solar  abundances 
to  illustrate  our  results  in  Figure  1 .  We  then  further  improve  the 
fit  by  varying  the  abundances.  We  find  that  in  order  to  better  fit 
the  sulfur  and  silicon  lines,  we  have  to  set  the  sulfur  to  20  times 
its  solar  abundance,  and  the  silicon  to  3  times,  while  all  the  other 
species  are  kept  at  solar  abundances.  The  x^  value  decreases 
from  0.3348  to  0.3 129  (by  about  7%).  While  some  of  the  lines 
are  better  fitted  (such  as  S  iv  1063  and  1073  A,  Si  iii  ^1010  A, 
and  Si  iv  ~1025  A),  the  S  iv  (1006  and  1 100  A)  are  far  too  deep 
(see  Fig.  2). 

4.2.  Accretion  Disk 

Next,  we  fit  the  solar  abundance  disk  models.  We  find  that  a 
low  inclination  is  needed  in  order  to  match  the  absorption  lines, 
and  in  our  models  we  initially  set ;  =  5°,  8°,  12°,  and  18°.  Again, 
because  we  have  no  information  about  the  mass  of  the  WD,  for 
each  value  of  log  g  (ranging  between  7.5  and  9.0)  we  vary  the 
mass  accretion  rate  between  10^'®  ^  and  10^*  °  Mq  yr^'  to  find 
the  best-fit  models.  Again,  we  find  that  the  best-fit  models  are 
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TABLE  2 

Synthetic  Spectra 


logs' 

(cgs) 

^WD 

(10^  K) 

Trot  sin  i 
(km  s“') 

[Si] 

(solar) 

[S] 

(solar) 

[Z] 

(solar) 

(deg) 

logM 
(Mo  yr-') 

WD/Disk 

(%) 

d 

(pc) 

xl 

Figure 

7.50 . 

34.0 

400 

1.0 

1.0 

1.0 

100/0 

377 

0.3427 

8.00 . 

36.0 

400 

1.0 

1.0 

1.0 

100/0 

266 

0.3361 

8.30 . 

37.0 

400 

1.0 

1.0 

1.0 

100/0 

247 

0.3348 

1 

8.30 . 

37.0 

400 

3.0 

20. 

1.0 

100/0 

211 

0.3129 

2 

8.50 . 

38.0 

400 

1.0 

1.0 

1.0 

100/0 

193 

0.3345 

8.50 . 

38.0 

400 

3.0 

20. 

1.0 

100/0 

184 

0.3100 

8.65 . 

38.0 

400 

1.0 

1.0 

1.0 

100/0 

172 

0.3360 

9.00 . 

40.0 

400 

1.0 

1.0 

1.0 

100/0 

135 

0.3390 

7.50 . 

1.0 

1.0 

1.0 

12 

-8.0 

0/100 

878 

0.3346 

7.88 . 

1.0 

1.0 

1.0 

12 

-8.5 

0/100 

878 

0.3346 

8.30 . 

1.0 

1.0 

1.0 

8 

-9.0 

0/100 

693 

0.3297 

8.30 . 

3.0 

20. 

1.0 

8 

-9.0 

0/100 

665 

0.3000 

3 

8.65 . 

1.0 

1.0 

1.0 

5 

-9.5 

0/100 

496 

0.3332 

9.00 . 

1.0 

1.0 

1.0 

5 

-10.0 

0/100 

361 

0.3376 

8.30 . 

32.0 

400 

3.0 

20. 

1.0 

08 

-9.0 

10/90 

700 

0.2989 

4 

8.30 . 

37.0 

400 

3.0 

20. 

1.0 

80 

-9.0 

74/26 

246 

0.2873 

5 

8.65 . 

38.0 

400 

3.0 

20. 

1.0 

80 

-9.5 

75/25 

184 

0.2880 

around  log  g  =  8.3,  but  the  ditFerence  between  the  xl  values  is 
not  significant,  and  the  improvement  over  the  single- WD  mod¬ 
els  is  also  only  marginal  (on  the  order  of  ~  1  %  in  the  xl  value). 
All  these  models  are  presented  in  Table  2.  Next,  we  improve 
the  disk  model  by  varying  the  abundances,  although  we  do  not 


expect  the  fit  to  improve  much  because  of  the  Keplerian  ve¬ 
locity  broadening  of  the  disk.  Flowever,  we  find  that  setting  the 
silicon  and  sulfur  abundances  to  3  and  20  times  solar,  respec¬ 
tively,  actually  reduces  the  xl  value  by  10%,  which  is  more  than 
for  the  single- WD  model.  This  is  because  the  S  iv  (1006  and 
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Fig.  3. — Nonsolar  abundance,  single-disk  model.  The  FUSE  spectrum  of  BB  Dor  is  shown  together  with  one  of  the  best-fit  synthetic  disk  models  {solid  black  line). 
The  model  hasM  =  0.80  Mq,  M  —  10“^  Mq  yr“\  i  —  8°,  a  distance  d  —  665  pc,  and  =  0.3000.  Abundances  and  solid  gray/dashed  gray  spectral  segments  are  as 
in  Fig.  2.  The  low  inclination  is  needed  in  order  to  fit  the  silicon,  sulfur,  and  carbon  absorption  lines.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this 
figure.] 
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Fig.  4. — Low-inclination,  nonsolar  abundance,  WD+disk  model.  The  FUSE  spectrum  of  BB  Dor  is  shown  together  with  one  of  the  best-fit  WD+disk  models  {solid 
black  line),  assuming  a  low  inclination  angle.  The  model  hasM  =  0.80  Mq,  M  —  10“^  Mq  yr“\  i  —  8°,  a  distance  d  —  700  pc,  and  =  0.2889.  Abundances  and  solid 
gray/dashed  gray  spectral  segments  are  as  in  Fig.  2.  The  WD  {dotted  line)  contributes  10%  of  the  flux,  while  the  high-inclination  disk  {dashed  black  line)  contributes  the 
remaining  90%.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 


1 100  A)  lines  have  better  fits  (Fig.  3).  This  best-fit  single-disk 
model  has  i=  8°,  M  =  10^^  Mq  a  distance  of  665  pc, 
and  xl  =  0.3000. 

4.3.  Composite  Model:  WD+ Accretion  Disk 

Last,  we  fit  composite  WD-i-Disk  models.  Since  the  number 
of  models  increases  exponentially  when  we  change  log  g,  Tej, 
and  M,  we  restrict  our  search  using  our  best  fit  for  log  g,  namely 
8.3.  We  also  set  the  sulfur  to  20  times  solar  and  silicon  3  times 
solar.  This  procedure  reduces  the  number  of  free  parameters  to 
3  (M,  Teff,  /). 

Low  inclination. —  Since  the  system  is  believed  to  have  a  low 
inclination  angle  (Chen  et  al.  2001),  we  generate  low-inclination 
models  (/  =  5°,  8°,  12°,  and  18°).  Not  surprisingly,  the  best  fit 
model  has  i  =  8°  and  M  =  10^^  Mq  yr^',  but  this  time  the  WD 
has  a  temperature  of  32,000  K  and  provides  only  10%  of  the 
flux,  while  the  disk  provides  the  remaining  90%.  Such  a  model 
again  brings  an  insigniflcant  improvement  in  the  value  of  and 
it  is  clear  (Fig.  4)  that  it  is  barely  distinguishable  from  the  best-fit 
single-disk  model.  From  the  point  of  view  of  the  physics,  this 
model  is  preferred  because,  if  the  system  has  a  low  inclination 
angle,  then  the  emission  ifom  the  WD  must  contribute  to  the  flux. 
Models  with  a  WD  temperature  T  <  30,000  K  have  a  WD  con¬ 
tribution  of  only  a  few  percent  of  the  total  FUV  flux,  and  cannot 
be  distinguished  from  the  single-disk  models.  If  the  WD  has  a 
temperature  7),^  <  30,000  K,  it  will  not  be  detected  while  the 
system  is  in  a  high  state  (with  M  =  10^^  Mq  yr^^).  These  results 
imply  that  the  contribution  from  the  WD  is  not  very  large  and 
that  the  temperature  of  the  WD  must  be  <32,000  K. 


High  inclination. —  From  the  sharp  Balmer  emission  lines, 
Chen  et  al.  (2001)  suggest  that  the  system  has  a  low  inclination. 
Flowever,  we  cannot  confirm  the  inclination  directly  from  the 
emission  lines  of  the  FUSE  spectrum,  as  all  the  sharp  emis¬ 
sion  lines  in  the  FUSE  spectrum  are  of  heliocoronal  origin.  Our 
low-inclination,  single-disk  models  provide  a  slightly  better  fit 
than  the  single-WD  models  (and  a  much  better  fit  than  the  high- 
inclination,  single-disk  models);  however,  for  completeness  we 
include  here  the  results  ifom  the  composite  WD-i-disk  model  fits 
when  the  assumption  about  the  inclination  angle  is  relaxed. 

We  search  for  the  best-fit  WD-l-disk  models  (assuming  log  g  = 
8.3)  in  the  T^s-M  parameter  space  using  all  inclination  angles, 
and  find  that  the  models  with  an  intermediate  inclination  (/  = 
18°,  41°,  60°,  and  75°)  do  not  provide  the  best  fit.  The  best-fit 
model  has  i  =  80°  and  reflects  a  situation  in  which  the  WD  is 
dominant  with  T  =  37, 000  K  and  contributes  3/4  of  the  total 
flux,  while  the  disk  has  M  =  IO^^Mq  yr^*  and  contributes  only 
1/4  of  the  flux.  The  distance  obtained  from  this  model  is  246  pc 
and  xj  =  0.2873.  This  is  the  least  x^  value  we  obtained  ifom  all 
our  models.  This  model  is  presented  in  Figure  5.  Similar  results 
were  obtained  assuming  log  g  =  8.65,  but  with  a  slightly  lower 
mass  accretion  rate  (Table  2). 

5.  DISCUSSION  AND  CONCLUSION 

BB  Dor  is  a  little-known  southern  NL,  and  consequently  both 
the  distance  and  the  mass  of  the  WD  are  unknown,  which  im¬ 
plies  a  larger  uncertainty  in  the  results.  In  addition,  the  FUSE 
spectrum  is  definitely  of  poor  quality.  In  theory,  a  fine-tuning  of 
the  temperature  (say  to  an  accuracy  of  about  ±50  K)  and  mass 
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Fig.  5. — High-inclination,  nonsolar  abundance,  WD+disk  model.  The  FUSE  spectrum  of  BB  Dor  is  shown  together  with  one  of  the  best-fit  WD+disk  models  {solid 
black  line).  The  model  hasM  =  0.80  Mq,M  =  10“^  Mq  yr“\  a  distance  ofd  —  246  pc,  andx^  =  0.2873.  The  inclination  angle  of  the  system  is  kept  as  a  tree  parameter, 
and  the  best-fit  model  gives  i  —  80°.  Abundances  and  solid  gray/dashed  gray  spectral  segments  are  as  in  Fig.  2.  The  fits  of  the  silicon,  sulfur,  and  carbon  absorption  lines  are 
obtained  by  letting  the  WD  component  dominate  the  spectrum.  The  WD  {dotted  line)  contributes  74%  of  the  flux,  while  the  high-inclination  disk  {dashed  black  line)  con¬ 
tributes  the  remaining  26%.  [See  the  electronic  edition  of  the  Journal  for  a  color  version  of  this  figure.] 


accretion  rate  can  be  carried  out  by  fitting  the  flux  levels,  such 
that  the  distance  to  the  system  (when  known)  is  obtained  accu¬ 
rately.  However,  the  fitting  to  the  distance  depends  strongly  on 
the  radius  (and  therefore  the  mass)  of  the  WD.  We  discuss  below 
some  additional  restrictions  that  we  use  to  constrain  the  proper¬ 
ties  of  the  system. 

On  the  basis  of  the  least  xj,  the  best  model  is  the  high-inclination, 
WD-Hdisk  composite.  However  (see  Fig.  5),  because  of  the  high 
inclination,  the  disk  contributes  a  rather  flat  component  all  the 
way  into  the  shorter  wavelengths.  At  such  a  high  inclination,  one 
would  not  expect  the  WD  to  dominate  the  flux,  but  rather  the 
WD  would  be  almost  completely  masked  by  the  swollen  disk; 
actually,  the  system  would  likely  be  observed  to  undergo  eclipses, 
but  none  have  been  observed.  Also,  since  BB  Dor  can  vary  as 
much  0.4  mag  in  1  hr,  this  means  that  the  light  cannot  be  dom¬ 
inated  by  the  WD,  and  that  the  disk  must  be  contributing  at  least 
40%  of  the  light.  In  other  words,  while  the  least  indicates 
that  the  best  fit  is  a  WD  with  a  rather  flat  disk  component,  there 
are  other  indications  that  this  cannot  be  correct.  Actually,  the 
WD-Hdisk,  high-inclination  model  strikingly  resembles  the  sec¬ 
ond  component  observed  in  the  FUSE  spectra  of  some  dwarf 
novae  during  quiescence  (e.g.,  VW  Hyi;  Godon  et  al.  2004).  This 
is  likely  a  caveat  in  state-of-the-art  spectral  modeling  rather  than 
an  indication  of  a  physical  link  between  the  spectrum  of  BB  Dor 
and  that  of  a  DN  in  quiescence.  The  need  for  improved  modeling 
also  stems  from  the  difficulty  in  producing  a  model  that  fits  low- 
and  high-order  ionization  lines  at  the  same  time.  For  BB  Dor,  it  is 


possible  that  the  S  iv  1063  and  1073  A  absorption  lines  form  in 
the  same  hotter  region /layer  where  the  O  vi  doublet  forms,  while 
the  C  III  and  Si  iii  lines  form  in  a  cooler  region/ layer  where  the 
Lyman  series  (and  continuum)  form.  This  is  similar  to  the  Hubble 
Space  Telescope  STIS  spectrum  of  TT  Crt  (Sion  et  al.  2008), 
which  exhibits  a  rich  variety  of  absorption  lines  from  different 
ionization  stages,  suggesting  line  formation  in  (at  least)  two  dif¬ 
ferent  temperature  regions. 

In  order  to  reduce  the  size  of  the  domain  for  which  we  have 
best-fit  models  in  the  parameter  space,  we  use  the  infrared 
magnitudes  J,  H,  and  K  from  the  Two  Micron  All  Sky  Survey 
(2MASS)  to  assess  the  distance  to  BB  Dor  as  prescribed  by  Knigge 
(2006)  for  systems  with  R  <  6  hr.  The  IR  data  were  collected  on 
1999  November  9,  at  a  time  when  BB  Dor  was  in  a  high  state 
with  a  visual  red  magnitude  R  =  14.60  and  a  blue  magnitude 
B=  13.90  (whereas  .5  ^  16.5  in  the  low  state).  The  2MASS  IR 
apparent  magnitudes  are  J  =  14.322,  H  =  14.089,  and  K  — 
14.053.  For  a  primary  star  with  A/wd  =  0.75  Mq  and  period  of 
3 .559  hr  (corresponding  to  BB  Dor),  the  donor  star  mass  is  M2  = 
0.25  Mq  (Patterson  et  al.  2005b  suggest  M2  =  0.256  Mq),  and 
the  IR  absolute  magnitude  estimates  are  A/j  =  1  Al ,Mh  =  6.90, 
and  Mk  =  6.63.  Inserting  these  into  equation  (15)  of  Knigge 
(2006)  gives  distances  of  235,  274,  and  305  pc,  respectively. 
These  distances  are  typically  underestimated  by  factors  of  2.05, 
1.86,  and  1.75  for  the  J,  //and  and  K  bands,  giving  distances  of 
482,  510,  and  534  pc,  respectively,  assuming  that  the  donor  star 
contributes  only  ~l/4-l/3  of  the  total  IR  flux  (Knigge  2006). 
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Fig.  6. — Dwarf  nova  systems  and  VY  Scl  NLs  above  the  period  gap.  The  sys¬ 
tems  are  plotted  in  the  (Teff,  P)  parameter  space  on  a  logarithmic  scale.  The  data 
for  the  NLs  are  taken  from  Hamilton  &  Sion  (2008)  and  references  therein.  The 
data  for  the  dwarf  nova  systems  are  taken  from  Sion  et  al.  (2008)  and  references 
therein.  For  RU  Peg,  we  have  used  the  more  recent  estimate  from  Godon  et  al. 
(2008).  The  data  point  for  V794  Aql  is  an  estimate  of  the  WD  temperature  based 
on  an  lUE  specttum  while  the  system  was  in  a  lower  state  (Godon  et  al.  2007). 
While  the  data  point  of  BB  Dor  is  lower  than  those  of  the  other  VY  Scl  systems, 
there  is  still  a  clear  separation  between  the  VY  Scl  systems  and  the  DN  systems, 
as  indicated  by  the  slanted  line.  [See  the  electronic  edition  of  the  Journal  for  a 
color  version  of  this  figure.^ 

The  distance  to  BB  Dor  is  therefore  certainly  larger  than  300  pc, 
and  most  likely  in  the  range  of  ~500  pc.  Since  BB  Dor  was 
observed  in  a  high  state  (with  2MASS),  it  is  likely  that  the  in¬ 
frared  flux  of  the  donor  star  is  on  the  lower  side,  namely  even 
less  than  1/4,  in  which  case  the  500  pc  itself  is  only  a  lower  limit. 
Using  this  restriction  on  the  distance,  we  can  safely  reject  the 
single- WD  models  (with  d  <  400  pc),  the  high-inclination  mod¬ 
els  (with  d  <  250  pc),  and  the  large-mass  log  g  —  9.0  models 
(with  d  =  361  pc  or  less). 

On  the  basis  of  the  distance  alone,  the  low-inclination  disk 
and  WD+disk  models  are  the  best  fits.  These  are  also  the  best 
fits  on  the  basis  of  the  optical  spectrum  (low  inclination)  and 
visual  magnitude  (high  state/disk). 

We  are  now  able  to  summarize  the  basic  characteristics  of  the 
system  and  its  WD.  We  are  confident  that  the  system  must  have  a 
WDwithlogg=  8.3  ±  0.3,7;ff  =  32,000 K or /ower(disk-tWD 
composite  models),  a  mass  accretion  rate  on  the  order  of  M  = 
10^^  Mq  yr^*  (disk  and  disk+WD  models),  and  a  distance 
somewhere  between  500  and  700  pc  (choosing  the  least  mod¬ 
els  agreeing  with  the  distance).  From  the  best-fit  single-disk 
models  alone,  we  find  that  the  inclination  angle  must  be  very 
small  (~8°)  to  fit  the  absorption  lines.  There  is  an  indication  that 
the  sulfur  and  silicon  are  overabundant,  and  it  is  also  possible 
that  some  of  the  high-order  ionization  lines  form  in  a  layer  in 
front  of  the  disk  and  the  WD. 


Although  the  FUSE  spectrum  of  BB  Dor  is  rather  poor,  and 
multiwavelength  data  are  limited,  we  have  narrowed  down  the 
region  in  the  parameter  space  {M,  P,  Teff )  for  this  system.  BB  Dor 
is  now  the  sixth  VY  Scl  NL  variable  with  a  temperature  esti¬ 
mate  for  its  WD.  With  a  temperature  of  <32,000  K,  the  WD  of 
BB  Dor  marks  the  lower  end  of  the  temperature  distribution  for 
WDs  of  VY  Scl  NL  variables,  which,  so  far,  for  the  other  five 
systems,  ranged  between  40,000  and  47,000  K  (Hamilton  & 
Sion  2008).  The  FUSE  spectrum  of  BB  Dor  clearly  shows  the 
drop  in  flux  in  the  shorter  wavelengths,  which  is  the  signature 
of  a  moderate  temperature  when  compared  to  the  FUSE  spec¬ 
tra  of  the  other  VY  Scl  systems  (see  e.g.,  V794  Aql;  Godon 
et  al.  2007).  Even  if  we  consider  the  best-fit  single- WD  model, 
the  WD  temperature  of  BB  Dor  only  reaches  37,000  K  (the 
model  with  40,000  K  and  log  g  =  9.0  gives  an  unrealistically 
close  distance  of  135  pc,  in  disagreement  with  the  expected 
IR  emission  from  the  secondary).  In  that  respect,  the  data  point 
in  the  (Tgjf,  P)  plane  for  BB  Dor  falls  under  the  region  of  the 
VY  Scl  NL  variables,  seemingly  into  the  region  of  the  param¬ 
eter  space  populated  with  dwarf  novae  (Sion  et  al.  2008).  How¬ 
ever,  as  shown  in  Figure  6,  the  data  point  for  BB  Dor  does  not 
seem  to  stand  apart,  and  the  separation  between  VY  Scl  systems 
and  DNs  can  be  made  easily  by  drawing  a  diagonal  line.  In  that 
respect,  it  is  actually  V794  Aql  that  divides  the  two  regions  of  the 
graph,  as  indicated  by  the  slanted  line. 
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